Experimental studies of the transient fluctuation theorem using hquid crystals 
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In a thermodynamical process, the dissipation or production of entropy can only be positive or 
zero according to the second law of thermodynamics. However the laws of thermodynamics are 
applicable to large systems in the thermodynamic limit. Recently a fluctuation theorem known as 
the Transient Fluctuation Theorem (TFT) which generalizes the second law of thermodynamics 
even for small systems has been proposed. This theorem has been tested in small systems such as 
a colloidal particle in an optical trap. We report for the first time an analogous experimental study 
of TFT in a spatially extended system using liquid crystals. 

PACS numbers: 05.70.Ln, 61.30.Gd 
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The laws of thermodynamics describe the physical be- 
haviour of macroscopic systems. The second lavif of ther- 
modynamics states that when such a system is taken 
from one equilibrium state to another in a process, the 
change in entropy can only be positive or zero depend- 
ing on wheather the process is irreversible or reversible 
respectively. Though the laws of thermodynamics are 
applicable under very general conditions, these laws are 
strictly valid only for large systems in the so-called ther- 
modynamic limit. For these large systems, the effects of 
thermal noise on the average macroscopic physical quan- 
tities are not manifested except under special physical 
conditions such as near phase transitions. However, when 
the system size is small or more precisely the change in 
the relevant energy of the system in a process is of the 
order of the thermal energy ksT, fc^ being the Boltz- 
mann constant and T being the absolute temperature of 
the system, the thermal noise is expected to play an im- 
portant role on it's behaviour. In particular, the validity 
of the second law of thermodynamics for small systems is 
of considerable debate since the time of Boltzmann. Re- 
cently a nonequilibrium fluctuation theorem (FT) known 
as the Transient Fluctuation Theorem (TFT) has been 
proposed to generalize the second law of thermodynam- 
ics for these small systems In its most general form, 
TFT not only predicts transient violation of second law 
of thermodynamics when the dissipation is comparable to 
the thermal energy fc^T but also it provides an expres- 
sion for the probability that a dissipative flux flows in 
a direction opposite to that required by the second law 
of thermodynamics. More precisely for a thermostated 
system at temperature T, this theorem states that in a 
time interval r, the probability P{U,t) of a dissipation 
being positive and the probability P{—^r) of same 
dissipation being negative in an irreversible process 
satisfies the following condition 
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The dissipation being an extensive quantity, the total dis- 
sipation increases as either the system size or the obser- 



vation time is increased. Then the above theorem implies 
that the production of entropy or positive dissipation will 
be overwhelmingly more likely than the consumption of 
entropy or negative dissipation in an irreversible process 
for large systems or for large observation time in accor- 
dance with the second law of thermodynamics. In this 
way the TFT generalizes the second law of thermody- 
namics even for small systems. Another FT for these 
small systems has also been proposed which relates the 
work done along nonequilibrium trajectories, in taking 
a system from one equilibrium state to another, to the 
thermodynamic free energy difference between the ini- 
tial and final equilibrium states of the systemQ- The 
physics of such small systems has recently gained wide 
interest in the scientific community Examples of such 
system include nano-materials, biological molecular ma- 
chines, quantum dots etc. In particular, the physics of 
such small systems when driven out of equilibrium is 
of paramount importance in many technological appli- 
cations of these systems. 

Though considerable theoretical and simulation 
studiesjj, d, 0, 0, S, 01 on these FTs have been reported 
in the literature, only a few experimental studies on the 
validity of these theorem have been performed. The main 
difficulty in the experimental studies of these theorems 
is in arranging the change in the relevant energy of the 
system under study to few fcsT required for the system 
to show deviation from macroscopic laws of thermody- 
namics. Nevertheless, the validity of these fluctuation 
theorems have been probed experimentally in small sys- 
tems such as colloidal particle in an optical trap 13, ll[ , 
in two- level systemH^, torsion pendulum immersed in a 
viscous fluid(13.] and deformation dynamics of single RNA 
molecule when stretched Q- In fact these theorems have 
been used to measure the folding free energy of a single 
RNA molecule when stretched [15] by an external applied 
force. However, all these experiments are performed on 
relatively small systems. The validity of these laws in 
the case of spatially extended systems where the physical 
properties are described by an effective order parameter 
which varies in space has not yet been probed experi- 
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mentally. In this letter, we report on the experimental 
studies of the TFT for such spatially extended systems 
using nematic liquid crystals. 

Liquid crystals (LC) usually made of rod like organic 
molecules lack the three-dimensional periodicity of crys- 
tals, but have anisotropic physical properties^). The 
simplest LC viz. Nematic LCs, which are used in practi- 
cally all commercial LC displays (LCD) have a long range 
orientational order of the long axes of the molecules. The 
orientational order described by the director n can be eas- 
ily deformed using external perturbations such as electric 
and magnetic fields. The distortions free energy of the 
director field can be described in terms of three elas- 
tic constants corresponding to the splay, bend and twist 
types of distortions of the director field[16| . In our exper- 
iments, a nematic liquid crystal is sandwiched between 
two appropriately treated tin oxide coated glass plates. 
The transparent tin oxide coating acts as electrode for 
the application of an electric field without obstructing the 
optical measurements. The glass plates are coated with 
polyimide and rubbed along a certain direction. These 
treatments of the glass plates align the long axes of the 
rod-like molecules parallel to the rubbing direction giv- 
ing rise to the so called homogeneous alignment of the 
liquid crystals between the plates. When such homoge- 
neously aligned nematic liquid crystals having positive 
dielectric anisotropy are subjected to an applied electric 
field perpendicular to glass plates, they exhibit a second 
order transition to a distorted structure (fig. [1^) above a 
threshold electric field. This is known as the Frederickz 
transition Above the Frederickz threshold, the equi- 
librium structure of the nematic liquid crystal depends 
on the applied electric field. 

As liquid crystals have anisotropic physical properties, 
the refractive index rig for light polarized parallel to the 
director n (extraordinary ray) is different from that of 
Uo corresponding to light polarized perpendicular to n 
(ordinary ray). This anisotropy in the refractive index 
of liquid crystals in this geometry gives rise to a phase 
difference ($) between the extraordinary and ordinary 
ray of the light beam passing through the sample and 
incident normal to the glass plates. The phase difference 
<i> depends on the orientations of the director n between 
the plates. In the absence of the applied electric field, the 
director n is parallel to the x-axis and $ = 27r And/A for 
light propagating along the z-axis, where An — rie ~ Uq 
is the birefringence of the LC sample and d, A are the 
thickness of the LC sample and the wavelength of the 
light respectively. When the applied electric field is in- 
creased beyond the Frederickz threshold, the director n 
progressively becomes parallel to the electric field and <& 
decreases with increasing electric field. Thus the change 
in the director configuration can be monitored by mea- 
suring the change in <I> using an interferometric technique 
under crossed polarizers. We have experimentally stud- 
ied the dissipative relaxation dynamics of $ when the 
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FIG. 1: Schematic representation of the experimental setup, 
(a) The distorted director configuration in nematic liquid 
crystal between two plates for applied electric field above 
Frederickz threshold. The dashed lines represent the local 
orientation of n which makes an angle 9 with respect to the 
X-axis. The double arrow represents the rubbing direction on 
the plates, (b) Schematic experimental setup where the He- 
Ne laser beam passes through a neutral density filter (NDF) 
and then through the LC sample held between crossed polar- 
izers (P) and (A). The intensity of the light passing through 
the sample is monitored using the photodiode (PD). 



liquid crystal is driven from one equilibrium state to an- 
other by a step change in the applied electric field. The 
change in the applied electric field is made sufhciently 
small such that the free energy difference between the 
initial and final state of the system is comparable to the 
thermal energy ksT and consequently study the validity 
of the TFT in this system. 

The schematic experimental setup is shown in fig. [T]d. 
The homogeneously aligned liquid crystal cell is placed 
between crossed polarizers with the director n making an 
angle 45 deg with respect to the polarizer. The intensity 
of a stabilized He-Ne laser (A = 632. 8nm) beam passing 
through the sample is monitored using a high gain low 
noise photodiode. The output of the photodiode is digi- 
tized using a 16-bit data acquisition board at the rate of 
1 kHz. We use a bipolar square wave electric field of fre- 
quency 1 kHz to drive the sample from one equilibrium 
state to another by changing the amplitude of the electric 
field. Use of AC electric field minimizes the effects of the 
charge impurities invariably present in the liquid crystals 
without affecting the director dynamics. The sample is 
thermostated in a heater with temperature stability of 
bmK. 

The experiments are performed on nematic liquid crys- 
tal (5CB) for sample thickness d = 27.0fim. The experi- 
ments are carried out at 30° C which is about 5°C below 
the transition temperature (35.1°C) from the isotropic to 
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FIG. 2: The average and fluctuating relaxation dynamics of 
X when the amplitude of the applied square wave voltage is 
changed from Vi = 1.490F to Vf = 1AS8V. Averaging is 
done over 3000 experiments. The average dynamics is well 
described by a single exponential decay to the flnal value with 
the relaxation time 362ms. 



the nematic phase. For a given temperature, we hold the 
sample with the amplitude of the applied square wave 
voltage at a certain value (Vi) for 2 seconds, then the 
amplitude of the applied voltage is changed in a step to 
a slightly lower value (Vf) and held at that value for 
subsequent 2 seconds. The normalized phase difference 
X = ($ — where i'i is the equilibrium value of <!> 

at the initial voltage Vi, is monitored as it relaxes from 
it's initial equilibrium value of zero to the final equilib- 
rium value X/ a.t voltage Vf. This cycle is repeated for 
3000 times to find the probability distributions of the dis- 
sipation over an ensemble of experiments. Fig. [5] shows 
the average relaxation dynamics of x averaged over 3000 
experiments and the fluctuating dynamics of x for one 
of the experiments when the amplitude of the applied 
voltage is changed from V = 1.490F to Vf = 1.488V". 
The average dynamics of x is described by a single ex- 
ponential decay from the initial value to the final value 
with a characteristic relaxation time ~ 362ms. Thus 
the dynamical evolution of x is analogous to the relax- 
ation dynamics of a particle in a harmonic potential in 
the overdamped regime in the presence of thermal noise 
described by Langevin equation. The effective strength 
(k^) of the harmonic potential associated with x for a- 
given applied voltage can be determined by measuring 
the thermal fluctuation of x and by the application of 
the equipartition theorem — ksT /a"^ where cr^ is the 
variance of x- For this purpose, we sample x in the initial 
and final equilibrium states of the system for 40 seconds 
and determine fc^ from the measured values of the vari- 
ance (T^ in these states. As Vi — Vf — 0.0021/ is small 



FIG. 3: Histograms of the dissipation Q.r obtained from an 
ensemble of 3000 experiments at different r (a) 5ms, (b) 80ms 
and (c) 125ms. The bin size A = 1.0. The solid lines are the 
fit to the normal distributions, (d) The corresponding plot 
of S{Q,Ti) versus the the dissipation firi associated with the 
i-th histogram bin for r = 80ms. The solid line is the best 
linear fit with slope 1.02 which agrees very well with the TFT 
predicted value of one. 



in our experiments, the measured values of ct^ and hence 
that of fc^ in the initial and final state of the system are 
found to be equal within the experimental error. The 
measured value of the variance = 1.13 X lO"'^ for 
1.490y at 30°C. 

The evolution of x after the step-like change in the 
applied voltage in this system is analogous to the over- 
damped Langevin type dynamics of a particle from the 
initial equilibrium position zero to final equilibrium posi- 
tion X/, when the harmonic potential felt by the particle 
is shifted from position zero to x/j keeping the strength 
of the potential fc^ same. Using the above analogy, the 



dissipation in time t can be written as 111 |17| 



knT 



iXr - Xo)- 
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where Xo and Xt are the values of x at time zero and t 
respectively. As expected, the expression for the dissipa- 
tion Qr is same as the work done in time r by a constant 
force fc^x/ on a particle which is in a harmonic poten- 
tial of strength fc,^ [1, @] . We calculate the dissipation fir 
from Eq. ([2]) for the 3000 repetitions of the experiment 
to get the desired histogram of ri-r- The histograms of 
ri-r for T — 5ms, 80ms and 125ms are shown in fig. [3l 
For r small, the average dissipation is small. Therefore 
the probability distribution of fir is peaked near zero 
with considerable probability of observing negative dis- 
sipation as can be seen from fig. [3^. As r increases, the 
average dissipation increases and the peak of the prob- 
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FIG. 4: The variation of the slope a with observation time 
T. Q is close to the TFT predicted value of one for r close 
to 80ms. The solid line is the best linear fit described by the 
equation a = 8.12 x 10"^r(ms) + 0.36. 



and for different applied voltages and we have observed 
similar variation of a with r. The details of these exper- 
imental results will be published elsewhere. 

In conclusion we have experimentally studied the va- 
lidity of the TFT in the spatially extended system us- 
ing liquid crystals. The TFT is studied for the dissi- 
pative dynamics of a macroscopic order parameter viz. 
the orientational order of the liquid crystals when driven 
from one equilibrium state to another. We construct the 
appropriate dissipation function for this system and ex- 
perimentally determine the histograms of the dissipation 
over an ensemble of experiments. The experimental re- 
sults are compared with the predictions of the TFT. We 
find agreement between the experimental observations 
and the predictions of TFT only for particular values 
of the observation time t. Simulation studies of spatially 
extended systems may shed further insight on to the va- 
lidity of TFT in these systems. We thank Dr. Abhishek 
Dhar for some helpful discussions and comments on the 
subject. 



ability distribution of f2,- shifts toward positive values 
of (fig. fig. [5};) and the probability of observing 
negative dissipation decreases. For very large r, the dis- 
sipation is much larger than ksT and the probability of 
observing negative dissipation becomes negligible as re- 
quired by second law of thermodynamics. 

To test the validity of the TFT in this system, we eval- 
uate the function S{Vlr) = ln[P{Vlr) / Pi'-^r)] from the 
experimentally measured probability distribution func- 
tion of D.r. Then according to the TFT, should 
be a linear function of f2r i-e. S{^r) = afJ^ with slope 
a = 1.0 for all r. We find from the experimentally mea- 
sured histogram of VIt the probability count {Ni) for dis- 
sipation between Vt^-i — A/2 and + A/2 correspond- 
ing to the i-th histogram bin Vlri = iA, A being the 
bin size. Then the experimental analog of the function 
S{^ri) = ln{Ni/N-i). We find that the experimentally 
determined S{^lTi) can be fitted well by a linear function 
of riri for all T. In fig. [3ji, we show the plot of S{riri) 
verses the dissipation for t — 80tos. The straight line 
in fig. shows the corresponding linear fit to the data 
with slope a = 1.02 agreeing very well with the prediction 
of TFT. Though the experimental results agree very well 
with the prediction of TFT for r close to 80ms, we find 
that in this system the value of the slope a is not equal to 
one for all r but it increases linearly with r. Fig. 0] shows 
the variation of the slope a with r which can be fitted 
with a linear equation a = 8.12 x lO^'^T(ms) + 0.36. We 
have repeated the experiments at different temperatures 
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